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ABSTRACT 

We demonstrate that the aperture mass as a measure for cosmic shear closely approximates (to better than w 5%) 
the scaled and shifted power spectrum of the projected mass density. This cosmological weak-lensing information 
can thus be used to directly infer the projected matter power spectrum with high accuracy. As an application, we 
show that aperture-mass observations can be used to constrain the cosmic density parameter and the power-spectrum 
amplitude. We show that, for a particular example, it should be possible to constrain iio to within « ±27%, and Og 
to within sa ±8% using weak-shear data on a square-shaped field of 8° side length. 



> 



1. Introduction 

A new measure for cosmic shear, the aperture mass (Map), was 
recently proposed by Schneider et al. (1998; hereafter S98). It 
was shown there that (mJ,) is related to the power spectrum P^ 
J^ of the projected density fluctuations filtered with a narrow func- 
ly-v tion. Here, we demonstrate in Sect. 2 that an accurate approxi- 
T— I mation to {M'L) can be constructed which directly yields P^- In 
C^ Essence, this makes the projected matter power spectrum a di- 
^^ fectly observable quantity, so that no deconvolution algorithms 
On need to be invoked. We then use our approximation in Sect. 3 to 
^^ infer the cosmic density parameter Q.q and the power-spectrum 
(— I normalisation Og from simulated data. We present our conclu- 
O sions in Sect. 4. 
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Aperture mass and approximations 



2.1. Effective convergence 



^ The two-point statistics of the gravitational lens properties of the 
k>( large-scale structure can be described, to high accuracy, in terms 
rN of the power spectrum of an equivalent single lens plane matter 

distribution (see, e.g.. Kaiser 1998; S98; and references therein), 

which is given by 
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where I is the Fourier conjugate to the angle 9, w and /a'(w) 
are the comoving radial and angular-diameter distance, respec- 
tively, a — (1 +z)^' is the scale factor, P^{k,w) is the density- 
perturbation power spectrum, and W{w) is the weight function 
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which depends on the probability distribution G{w) of source 
distances. The upper integration limit wh is the comoving hori- 
zon distance, here defined as the comoving distance correspond- 
ing to redshift infinity. Quite intuitively, eq. (nl) relates the lens- 
ing power on angular scales 9 — 2kI^^ to the power in density 
fluctuations at a comoving scale Ink^^ ~ fK{w)Q. 

We parameterise the source-distance distribution G{w) as a 
function of redshift, G,;(z), specified by 
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It is normalised to < z < oo and provides a good fit to the ob- 
served redshift distribution (e.g. Smail et al. 1995). The mean 
redshift (z) is proportional to zq. For p = 1.5 which we assume 
throughout, (z) w 1.505zo- 

2.2. Aperture mass 

The aperture mass Map (9) as a function of smoothing scale 9 is 
defined in terms of a weighted average within a circle of radius 9 
of the surface mass density of the equivalent single lens plane. It 
can readily be obtained from the observed image ellipticities of 
faint background galaxies which provide an unbiased estimate 
of the shear y. 
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by replacing the integral over the shear by a sum over galaxy 
ellipticities. Here, Yt is the tangential component of the shear 
relative to the aperture centre, and Q is an appropriately cho- 
sen weight function which is non-zero only for < i3^ < 9- 
The mean-squared aperture mass is related to the effective- 
convergence power spectrum through 
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where 7(r|) is related to the filter function Q (see S98). For 
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and J4(r|) is the fourth-order Bessel function of the first kind. 
7^(r|) peaks at tjq « 4.11. Examples for the nns aperture mass 
(Map) '' ^ are plotted in Fig. g as a function of aperture radius 9. 
It is crucial in eqs. (|l]) and (^ to take the non-linear evolu- 
tion of the density power spectrum into account. For aperture 
radii of order a few arc minutes, 9 « 3 x lO^'^rad, the peak 
TJQ of J^{r\) translates to / « 1.4 x lO'*, which corresponds to 
Ink^^ « 1/i^^Mpc for sources around redshift unity, i.e. the 
physical scales of density perturbations to which the aperture 
mass is most sensitive to are well in the non-linear regime of 
evolution. We assume a CDM power spectrum and describe its 
non-linear evolution as given by Peacock & Dodds (1996). We 
normalise the spectrum to the local abundance of rich clusters 
by choosing Og as derived by Viana & Liddle (1996) and Eke, 
Cole &Frenk (1996). 



2.3. Signal-to-noise ratio 

Although an nns aperture-mass ampHtude around one per cent 
appears low, the signal-to-noise ratio of the aperture mass can be 
quite high. There are three sources of noise in a measurement 
of (Ml,). Since it will be inferred from distortions of galaxy 
ellipticities, the intrinsic non-vanishing ellipticity of the galaxies 
provides one source of noise, and the random positions of the 
galaxies provides another. The third source of noise is due to 
cosmic variance. Assuming a large number of galaxies A^ per 
aperture, and neglecting the kurtosis of the aperture mass, the 
dispersion of Map in a single aperture is (see S98) 
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where o^ ~ 0.2 is the intrinsic dispersion of the galaxy elliptic- 
ities. On angular scales of a few arc minutes and smaller, the 
galaxies dominate the noise, while the cosmic variance domi- 
nates on larger scales. Of course. Map will be measured in a 
large number of apertures A^ap rather than a single one. If the 
apertures are independent, the dispersion O) is reduced by a fac- 
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tor of N„o , and the ensemble variance becomes 
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Typical signal-to-noise ratios reach values of > 5 for aperture 
radii of a few arc minutes and data fields of a square degree in 
size. 

2.4. Approximate rms aperture mass 

The strong peak of 7^(r|) - see Fig. 2 in S98 - motivates the 
approximation of J^ by a delta function. 
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with A =512/(1 15571^) « 1.43 X lO^^, rjo = 6937i:/512«4.25, 
as determined from the norm and mean of 7^]^ Then the expres- 
sion for the aperture mass (Bl) becomes particularly simple, 

(M2p)(e) « {M%){Q) EE 6(57r)-ie-2pK(ilo/e) . (10) 

If this provides a good approximation, the observable rms aper- 
ture mass at angular scale 9 would directly yield the effective- 
convergence power spectrum Pk{1) at wave number r|o/9, and 
thus a most direct and straightforward measure for dark-matter 
fluctuations. We therefore have to investigate whether replacing 
the filter function J^ (r|) by a delta function can be justified. 

Let us first assume that P-^il) can be approximated locally as 
a power law in / over a range in which J^ differs significantly 
from zero, 

P^{l)=Bl"'=ff . (11) 

Then, the true aperture mass (pb becomes 
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'in fact, the peak of /-(t|) can well be approximated by a Gaussian. A least- 
square fit of a Gaussian to fi{r\) yields 

y2(ii) « A'exp[-(Ti-Tii)V(2a2)] 

with mean tIq « 4.11, amplitude A' a; 4.52 X 10"^, and variance a ks 1.24. The 
small width of the Gaussian encourages the approximation (bI). 



where the Gamma functions arise from integrating the power- 
law spectrum times the squared fourth-order Bessel function. 



The approximate aperture mass (M^p) defined in eq. (10) be- 
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The relative deviation of the approximation from the true aper- 
ture mass is a function of the effective power-law exponent n^ff 
only. It is plotted in Fig. H. Clearly, the deviation of (M^p) 
from (Map) is very small, less than five per cent, for effective 
power-spectrum slopes in the range — 1.5 < n^ff < 0.5, and this 
is exactly the range of slopes in the I interval contributing most 
of the power to the aperture mass for aperture radii of > 1'. It is 
therefore fair to say that, especially in the presence of measure- 
ment errors, the aperture mass and its approximation ( p^ can be 
considered equivalent. Figure ||shows (Map) and (M^p) for three 
cluster-normalised CDM models as examples for realistic non- 
power-law spectra, and emphasises the conclusions from Fig. nl 
The approximation ( [l(3| ) to the aperture mass is excellent. This 
immediately implies that measurements of the aperture mass di- 
rectly measure the effective-convergence power spectrum, and 
therefore the latter effectively becomes an observable quantity, 



p.{l)^{5K/6){^o/lf{Ml„){^o/l) 



(14) 



This provides the most straightforward access to a quantity of 
paramount importance for cosmology, i.e. the dark-matter power 
spectrum. 




Fig. 1. — Relative deviation between the true aperture mass {M^„) and 
the approximate aperture mass (mL) ^^ defined in eq. (UQ), for an 
assumed effective-convergence power spectrum that is a power law 
in /. For the most relevant range of effective power-law exponents, 
— 1.5 < «eff lS 0-5, the relative deviation is less than five per cent. 



3. Applications 

It is now interesting to investigate a specific application of our 
results. Let us assume we are given a large, quadratic data field 
of angular size L, into which we place apertures of radius 9. 
We can place apertures of a fixed size densely on the field since 
the correlation of the aperture mass in neighbouring apertures 
is very small due to the narrowness of the filter function J^{r\) 
(see Fig. 8 in S98). However, we want to sample (Map) at a 
variety of aperture sizes rather than a single one. Starting from a 
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Fig. 2. — Exact and approximate rms aperture masses (M^p) and (Map) 
are shown for three cosmological models. All are cluster-normalised 
CDM models in which the non-linear growth of the density power spec- 
trum ^5 was taken into account. Solid curves: low-density, open model 
(fl() = 0.3, flyv = 0, h = 0.7); dotted curves: low-density, spatially flat 
model (fio = 0.3, 0.\ = 0.7, h = 0.7); dashed curves: Einstein-de Sit- 
ter model (flo = I, Q.A = 0, h = 0.5). The redshift distribution (|) of 
the sources was used. Typical rms aperture masses reach --^ 1 per cent 
at angular scales of a few arc minutes. The normalisation to the local 
abundance of rich clusters causes the curves to be very similar for the 
different choices of the cosmological parameters. While the approx- 
imate aperture mass slightly underestimates the true aperture mass at 
small aperture radii, 6 > 1', by at most ~ 5%, the approximations fall 
exactly on top of the true curves for angular scales beyond ~ 10'. 



smallest aperture radius Gq, we choose the next largest aperture 

radius as \/29o and so forth, so that consecutive aperture radii 
encompass on average twice as many galaxies. Clearly, the data 
for different aperture sizes are not completely uncorrelated, but 
decorrelate quickly with differing radius and centre separation - 
see again Fig. 8 of S98. Therefore, we assume that the number 
of independent apertures of radius 9 is given by A^ap = [L/(29)]^, 
and the expected number of galaxies per aperture is A^ = nnQ^. 
We assume a galaxy density of n = 30 per square arc minute. 
The expectation value of (M^ )(9) is given by eq. (js]). Its l-o 
error at fixed aperture radius is expected to be given by eq. ^. 
We show examples of simulated measurements M, together with 
their expected l-o error bars in Fig. |] for field sizes of 4° and 
8°. 

Let Mi be the measurement of (M^p) for aperture radius 9,-, 
then (Q) yields the expected error of M,-. We can then try and 
fit the ensemble of measurements with the approximate aperture 
mass (Map). For that purpose, we define the x^ function 
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As an example, we assume that data points like those shown in 
Fig. Ijhave been measured in a cluster-normalised, low-density 
open CDM model. Assuming the CDM shape of the density- 
perturbation power spectrum, we can then vary parameters such 
as to minimise (llSb. The sensitivity of the aperture mass to 
the Hubble constant is very small. The remaining parameters 
are then the parameters p and zo of the source redshift distribu- 
tion, the amplitude of the power spectrum, parameterised by Og, 
the density parameter iio, and the cosmological constant Q.A- 
lt can safely be assumed that the redshift distribution G{w) of 



Fig. 3. — Simulated aperture-mass measurements employing the sam- 
pling strategy detailed in the text. The error bars show l-o errors 
derived from eq. M). Field sizes of 4° and 8° were assumed for the 
open squares with dotted error bars and the filled squares with solid 
error bars, respectively. The true aperture mass is shown as the solid 
curve. The low-density, open CDM model was assumed for this plot, 
and zo = 1 was chosen for the redshift distribution of the sources. Ac- 
cording to eq. (lid), the effective-convergence power spectrum can di- 
rectly be read offfrom this plot. 



background galaxies is known with high accuracy from other 
observations. We can therefore restrict ourselves to essentially 
two free parameters, ag and iio, and choose either Q.^ = or 
Q.A = 1 — iio- X^ contours for the former case are shown in 
Fig.|. 

The solid contours in the left panel of Fig. ^ show that at the 
l-o significance level, iio and Og can be constrained to relative 
accuracies of « ±27% and w ±8%, respectively. Table |l] gives 
l-o limits for fields of 4° and 8° side length. Similar experi- 
ments using Q.A = 1 — iio instead of Q,a = yield a minimum 
X~ > 10, so that spatially flat cosmologies can be rejected at very 
high significance using the data simulated with Q.a ~ 0. 

Of course, the quality of the results depends on how accu- 
rately we guess the shape of the true power spectrum. For the 
contours in Fig. 0, we vary two parameters and keep all others 
fixed. In Fig. ^for example, we keep /i = 0.7 constant, so that the 
shape parameter F = Q.oh of the power spectrum is F = 0.7iio. 
If we would set the shape parameter instead to F' = /F with 
/ ^ 1, the best-fitting parameter combination (iio, Og) would 
shift towards /iio along the "valley" indicated by the %^ con- 
tours without substantially degrading the quality of the fit. This 
degeneracy reflects the degeneracy in the cluster abundance used 
to normalise the power spectrum for simulating the input data, 
because (M^p)''^ essentially measures the halo abundance at in- 
termediate redshifts, z ~ 0.3. To illustrate this point, we indicate 
the present-day cluster-abundance constraint as the long-dashed 
curve in Fig. Q. It is seen to follow well the "valley" in the x^ 
contours for low Q.q, but it departs increasingly for increasing 
Q.0 because structure grows more rapidly between z ^ 0.3 and 
today for higher iio. The parameter degeneracy with the cos- 
mological constant can be broken using the skewness of Map, 
for which the degeneracy between iio and Q.^ is different (Van 
Waerbeke, Bernardeau & Melher 1999). 

Although our assumptions concerning the number of inde- 
pendent apertures that can be placed inside a given data field 
and the neglect of the kurtosis in the approximation for o(Ma„) 
may be slightly optimistic, we believe that these effects do not 
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In contrast to other methods for constraining Pk{1) through 
cosmic-shear data (Kaiser 1998; Seljak 1998; Van Waerbeke 
et al. 1999), the aperture-mass method allows for a simple 
combined analysis of independent data fields scattered across 
the sky. Furthermore, since (M^p)''^ is a scalar quantity di- 
rectly obtained from observed galaxy-image ellipticities, its full 
probability-distribution function can be derived and used for pa- 
rameter extraction. 

The approximate aperture mass can then straightforwardly 
be applied to constrain cosmological parameters from observed 
aperture masses. We simulated observations of (M^) and their 
expected errors and showed that Q.q and Og can be recovered 
with relative accuracies of w ±27% and w ±8%, respectively, 
using weak-shear data on a square field with 8° side length. For 
that, we have assumed that the parameters of the source redshift 
distribution are sufficiently well known. We believe that this is 
not a serious limitation because of the high reliability and accu- 
racy of photometric redshift determinations (e.g., Benitez 1998, 
and references therein). In fact, the source-redshift dependence 
of (Mg ) can be used as a consistency check: With increasing 

source redshift, (Mi,) increases, while the skewness of Map de- 
creases. 



Fig. 4. — Contours of x as defined in eq. ( |l5[ ) are shown in the do — Og 
plane. Input data and their errors were simulated as detailed in the 
text, assuming quadratic data fields of side length 4° and 8° (dotted 
and solid contours, respectively). Contours are shown for 1, 2, and 
3-a significance levels. The input cosmological model was the cluster- 
normalised, low-density, open CDM model with flo = 0.3, Q.\ = and 
h = 0.7. Due to the cluster normalisation, CTg = 0.85. 0.\ = was 
assumed for the fit. The contours show that data on a 8° x 8° field are 
sufficient to constrain Hq at the l-o level to within 0.18 < fiy < 0.34, 
and the amplitude to 0.85 < csg < 0.99. The long-dashed curve shows 
the present-day cluster-abundance constraint on Hq and Og . 



change our conclusions appreciably. It should be noted that for 
applying our technique, the data field need not be a connected 
region, l)ut arbitrarily placed smaller areas on the sky with about 
equal total solid angle are equally useful. The rapid evolution 
towards square-degree CCD cameras at the best observing sites 
leaves us highly optimistic concerning the future application of 
our method. 
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Table 1. — l-o Hmits to iio and og obtained with fields of 4° 
and 8° side length. 
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input 


fit 


input 


Og 

fit 


r 

8° 


0.3 
0.3 


0.25 ±0.12 
0.26 ±0.08 


0.85 
0.85 


0.93±0.12 
0.93 ±0.08 



4. Conclusions 



We have shown here that a recently proposed measure for cos- 
mic shear, the dispersion of the aperture mass (Map(9)), can be 
approximated by the local value of the power spectrum f k of the 
projected matter fluctuations at wave number r|o/6 to a relative 
accuracy better than 5%. This is possible because the aperture 
mass is a convolution of P^ with a narrow filter function which 
for this purpose can safely be approximated by a Dirac delta 
function. 



